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Abstract 

Measurements of from event shapes in e+e^ annihilation are discussed including 
recent determinations using experimentally optimized scales, studies of theoreti- 
cally motivated scale setting prescriptions, and recently observed problems with 
predictions in Next to Leading Logarithmic Approximation. Other recent precision 
measurements of ag are briefly discussed. The relevance of power terms for the 
energy evolution of event shape means and distributions is demonstrated. Finally 
a summary on the current results on as(M^) and its running is given. 
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1 Introduction 



Most of our current knowledge about the strong coupling and its running comes from 
the analysis of the process e"'"e~ hadrons which is the simplest possible initial state 
for strong interaction. The energy scale of the process is exactly known and hadronic 
interaction is limited to the final state. For many observables the hadronization corrections 
are proportional 1/Q and therefore LEP with an energy range up to about 200 GeV in the 
year 2000 is an ideal laboratory for studying QCD. Due to the huge cross section at the 
Z° resonance every LEP experiment collected several million hadronic events, serving for 
a precise determination of as(M|). With the high energy data from LEP 2, the running 
of a<j and the influence of non-perturbative power terms on the observed quantities can 
be studied in detail. 



2 Consistent Measurements of cXs from the Analysis 
of Data using Oriented Event Shapes 

Within the recent analysis of Z° data by the DELPHI collaboration, the distributions 
of 18 different shape observables are determined as a function of the polar angle i^T of the 
thrust axis with respect to the e"'"e~ beam direction. Since the definition of the thrust axis 
has a forward-backward ambiguity, cos?? > has been chosen, cos'd is called the event 
orientation. The definition of the observables studied can be found in [|^ . The theoretical 
predictions in C(o;^) are calculated using EVENT2 0, which applies the matrix elements 
of the Leiden Group [0]. Using this program, the double differential cross section for 
any IR- and collinear safe observable Y in e"'"e~ annihilation in dependence on the event 
orientation can be calculated: 



^ ' = asifi ) ■ A{Y, cos ^) (1) 



atot dYd cos 

+ ■ \ B{Y, cosd) + (27r/3oln(x^) -2)A(F,cos^9) 

where ols = as/^n and Po = (33 — 2nf)/127r. atot is the one loop corrected cross section 
for the process e"'"e~ — > hadrons. A and B denote the 0{as) and O^a"^) QCD coefficient 
functions, respectively. 

„2^ 



The dependence on the renormalization scale enters logarithmically in 0{a2)- The 
is defined by = a;^( 



scale factor is defined by fi"^ = a;^Q^ where Q = Mz is the center of mass energy. In 



0{a^), the running of the strong coupling at the renormalization scale fi is given by 

Pi Inlng" 



(5) 

where A = A^^ is the QCD scale parameter computed in the MS scheme for Uf = 5 
flavors and Pi = (153 — 19?7,j)/247r^. The renormalization scale /i is a formally unphysical 
parameter and should not enter at all into an exact infinite order calculation. However, 
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Figure 1: (a) QCD fits to the measured thrust distribution for two bins in cos-i^T- (b) 
Measured thrust distribution at various fixed values of 1 — T as a function of cos-i^x- The 
sohd hues represent the QCD fit. 

within the context of a truncated finite order perturbative expansion for any particular 
process under consideration, the definition of n depends on the renormalization scheme 
employed, and its value is in principle completely arbitrary. 

The traditional experimental approach is, to measure all observables at the same, fixed 
scale value, the so-called physical scale (PHS) x^ = 1 or equivalently fi"^ = Q^. Apply- 
ing PHS to the high precision DELPHI data at ^/s = Mz yields in general only a poor 
description of the measured event shape distributions, values up to I'^df ~ 40 are 
found. For the PHS choice the 2^^*^ order contribution in Eq. |l] is in general quite large. 
For some observables the ratio of the 0{a1) with respect to the 0{as) contribution r^vLo 
is almost Ir^Lol — 1-0, indicating a poor convergence behavior of the corresponding per- 
turbative expression. This quite naturally explains the observation of the wide spread of 
the measured as values for the individual observables (see Fig. ^d). If PHS is applied, an 
unweighted average for the as{M'^) values of the 18 observables yields x^ /ndf = 40/16, i.e. 
the individual measurements are inconsistent. For the differential 2-jet rate determined 
using the Geneva-Algorithm, the fit applying = 1 fails completely to describe the data. 
Here, no a, (M|) value can be derived at all if PHS is applied. Therefore, the central 
method for measuring a;s(M|) has been chosen to be a combined fit of as and the scale 
parameter x^, a method known as experimentally optimized scales (EXP). Here one finds 
in general much smaller contributions from the 0{a'^) term in Eq. |I], indicating a better 
convergence behavior of the perturbative series. Applying EXP, the O(a^) predictions 
including the event orientation yield an excellent description of the high statistics data 
(see Fig. |l| as an example). For all observables considered, the QCD fit yields x^ /ndf — 1. 

Fig. § shows the renormalization scale dependence of as{M'^) for some of the observ- 



3 



;:r 0.15 
S 

0.14 
0.13 
0.12 
0.11 



200 
150 
100 
50 



"-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 

Figure 2: as(M|) and Ax^ = ~ xLin from 0{al) fits for some of the observables 
studied as a function of the renormahzation scale x^. Additionally, the minima are 
indicated in the as(M|) curves. 

ables studied. It turns out that in order to describe the data, the scale has to be fixed to 
a rather narrow range of values. Consistent as(M|) measurements can only be derived, 
if the optimized scale values are applied, i.e. from the values corresponding to the 
minima of the individual fits. For most of the observables the scale dependence in 
the vicinity of the minima is significantly reduced, but even for the few observables 
exhibiting a strong scale dependence around the minima, the corresponding q;s(M|) 
values are perfectly consistent. The observable with the smallest scale dependence of 
a^j(M|) is the Jet Cone Energy Fraction (JCEF) [^]. Here, the change in as(M|) is only 
Aas{M^) = ±0.0010, even if the scale is varied within the large range of < < 1. 
Additionally JCEF has the smallest hadronization correction uncertainties as well. 

The as(M|) values determined from 18 different observables are shown in Fig. ^ 
for EXP in comparison with PHS. For EXP the scatter among the different observables 
is significantly reduced. The errors of as(M|) correspond to the quadratic sum of the 
uncertainty from the fit, the systematic experimental uncertainty and the hadronization 
uncertainty. Conservatively, an additional uncertainty due to the variation of the central 
renormahzation scale value between 0.5 • x^ and 2. ■ x^ has been considered for both 
methods PHS and EXP. An unweighted average yields a<,(M|) = 0.1165±0.0026 for EXP 
to be compared with a,(M|) = 0.1243 ± 0.0080 in the case of PHS. The corresponding 

value is x'^/^df = 7.3/17 for EXP and x'^/^df = 79/16 for PHS. It should be empha- 
sized, that the consistency for EXP does not depend on the additional uncertainty due 
to renormahzation scale variation. Ignoring this uncertainty yields an consistent average 
value of as(M|) as well {x'^/n^f = 9.2/17). A weighted average of q;s(M|) considering 
correlations between the observables yields as(M|) = 0.1164 ± 0.0025, almost identical 
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Figure 3: Results of the as(M|) measurements from 18 event shape distributions, (a) 
fits using experimentally optimized scale values, (b) fixed scale fits : = 1. The errors 
on as{Mz) indicated correspond to the quadratic sum of the uncertainty from the fit, the 
systematic experimental uncertainty and the hadronization uncertainty. The dotted error 
bars indicate the additional uncertainty due to the variation of the central renormalization 
scale value between 0.5 ■ and 2. ■ x^. 

to the unweighted average. The investigation of the influence of heavy quark mass effects 
on as(M|) is under study. A preliminary estimate yields as(M|) = 0.117 ± 0.003. 

2.1 Theoretically Motivated Scale Setting Methods 

Additional studies |jl| have been made using theoretically motivated scale setting pre- 
scriptions, like the principle of minimal sensitivity (PMS), the method of effective charges 
(ECH) and the method of Brodsky, Lepage and MacKenzie (BLM). In the case of PMS 
and ECH, a strong correlation with the measured scale values from EXP can be observed. 
For the BLM method no such correlation is observed (see Fig. Moreover, the BLM 
fits do not converge for some of the observables under study. The individual as values 
from the remaining observables turn out to be inconsistent. However, the average val- 
ues of as(M|) for all methods considered are consistent with EXP, but the scatter of 
the individual as(M|) measurements is somewhat larger for the theoretically motivated 
methods. 

The influence of higher order contributions has also been investigated |jl| by using the 
method of Pade Approximants for the estimate of the uncalculated 0{a^) contributions 
(PAP). In comparison with (9(a^), the renormalization scale dependence for PAP is sig- 
nificantly reduced. Here, a fixed scale value of = 1 has been chosen for the as(M|) 
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Figure 4: Correlation of experimentally optimized renormalization scale values with values 
predicted by theoretically motivated scale setting prescriptions, {left side: ECH method, 
right side: BLM approach) 

measurements from the individual observables. The average value of as{M'^) is again in 
perfect agreement with the 0{a1) result, suggesting small contributions due to missing 
higher order predictions. 

2.2 Comparison with NLLA predictions 

The probably most relevant check on the influence of higher order contributions comes 
from a study of the all orders resummed next to leading logarithmic approximation 
(NLLA), which has been calculated for a limited number of observables. Two differ- 
ent strategies have been applied Pure NLLA calculations have been used to measure 
as(M|) in a limited kinematical region close to the infrared limit, where the logarithmic 
contributions become large. Matched NLLA + 0{a1) calculations have been used to ex- 
tend the 0{al) fit range towards the 2-jet region. Unlike 0{a1) theory, no optimization 
is involved in adjusting the renormalization scale for NLLA predictions. Therefore the 
renormalization scale value has been fixed to = 1. Both methods yield average values 
of a;s(M|) compatible with the average from 0{a1), the scatter of the individual mea- 
surements is somewhat larger for the resummed theory. 

Looking at the individual measurements, one finds that the values derived from 
matched predictions are systematically higher than those from pure NLLA theory. For 
some of the observables, the as value from matched predictions are even higher than for 
both pure NLLA and 0{a1) predictions. Clearly the matched result is expected to be 
a kind of average of the results from both the distinct theories. Moreover, the matched 
theory predictions for some of the observables yield only a poor description of the high 
precision data, a more detailed investigation reveals a systematic deviation of the predicted 
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Figure 5: Comparison of the measured C-Parameter Distribution with three different 
QCD Fits: z) 0{al) using an experimentally optimized renormalization scale, ii) C(a^) 
using a fixed renormalization scale x^ = 1 and in) InR matched NLLA (x^ = 1) for the 
C Parameter. The lower part shows the relative difference (Fit-Data)/Fit. Whereas the 
C(a^) curve describes the data over the whole fit range, the slope of the curves for the 
fixed scale and InR matching show a similar systematic distortion with respect to the 
data. 

slope with respect to the data (see Fig. The distortion observed is similar to the 
distortion obtained using O(a^) predictions applying a fixed renormalization scale value 
of = 1, indicating that the terms included from the 2""^ order predictions dominate the 
matched predictions. Whereas = 1 seems to be an appropriate choice for pure NLLA 
predictions, the similarity of the two fit curves indicate a mismatch of the renormalization 
scale values for the 0{af) and NLLA part of the combined prediction. 

3 Other High Precision cXg Measurements 



3.1 cxg Determination from Precision Electroweak Measurements 

The precise electroweak measurements performed at LEP and SLD can be used to check 
the validity of the Standard Model and to infer information about its fundamental param- 
eters. Within the Standard Model fits of the LEP Electroweak Working Group, the value 
of a<i(M|) depends mainly on Ri, Tz and atad- The theoretical prediction is known in 
NNLO-Precision. A recent update of the fit results presented at the ICHEP'98 |^ yields 
a value for the strong coupling of as(M|) = 0.119 ± O.OOSfu ± 0.002Theo- 
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Figure 6: The vector and axial- vector spectral functions in hadronic rdecays. The data 
points correspond to the sum of all contributing channels. Some exclusive contributions 
are shown as shaded areas. 

3.2 as Determination from Spectral Functions in Hadronic r 
Decays 

One of the few processes, where QCD predictions are known in NNLO, is the hadronic 
decay of the r Lepton. Within the recent OPAL analysis , has been determined from 
the moments of the spectral functions of the vector and axial-vector current in hadronic r 
decays, which are weighted integrals over the differential decay rate dRry/A/ds for vector 
(V) and axial-vector (A) decays: 

Rr,V/AiSo) = jds(^l-y^ (^-^^ dRry/Alds (3) 


The analysis involves a measurement of the invariant mass of the hadronic system and 
thus requires an exclusive reconstruction of all hadronic final states. After unfolding the 
measured spectra, normalizing them to their branching ratios and summing them up with 
their appropriate weights, one obtains the spectral functions shown in Fig. |^. 

The QCD predictions have been calculated including perturbative and non-perturba- 
tive contributions. Within the framework of the Operator Product Expansion (OPE) the 
non-perturbative contributions are expressed as a power series in terms of In con- 

trast to the perturbative part, the power corrections differ for the vector and axial- vector 
part, thus the difference of the vector and axial- vector spectral function is sensitive to 
non-perturbative effects only. The perturbative contribution is known in ) and partly 
known in 0{a\). On a first glance this accuracy looks quite impressive, it should how- 
ever be emphasized, that with the small momentum transfers involved, the divergence 
of the perturbative series is a much more important issue for the r decay than for ag 
measurements at high energies. Recent estimates of the complete 0{af) contribution]^ 
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Prediction as(M|) xV^d/(Fit 1) xV^d/(Fit 2) 

FOPT 0.1191 ±0.0008exp±0.0013t;,eo±0.0003e™i 0.17/1 0.62/4 
CIPT 0.1219 ±0.0010exp±0.0017t;,eo±0.0003e™i 0.16/1 0.63/4 
RCPT 0.1169 ±0.0007exp±0.0015t;,eo±0.0003e™i 0.07/1 0.61/4 



Table 1: Results on as(M|) from an OPAL analysis of spectral functions in hadronic 
T decay. Given are the as{M'^) values for the three different perturbative calculations 
studied. The as{M^) values for the two different fits are identical. Shown are also the 
X^/^df values for both fit strategies. 

to the perturbative prediction by the means of Pade Approximation predict an O(a^) 
contribution of about 20 % w.r.t. the leading order term, indicating that further higher- 
order terms could have a significant effect on the perturbative prediction. Three different 
calculations for the perturbative prediction have been studied within the OPAL analysis. 
Apart from the standard fixed order perturbative expansion (FOPT), two attempts have 
been studied to obtain a resummation of some of the higher order terms. The so-called 
contour improved perturbative theory (CIPT) accounts on higher order logarithmic terms 
in as by expressing the perturbative prediction by contour-integrals, which are evaluated 
numerically using the solution of the renormalization group equation (RGE) to four-loops. 
The third calculation includes an all order resummation of renormalon contributions in 
the so-called large /3o-hmit (RCPT). This strategy has the advantage to be renormaliza- 
tion scheme invariant. 

Two different fits to the data have been performed. The first fit uses 5 moments 
from the sum of the vector and axial-vector spectral function for the determination of 
as(m^) together with three parameters from OPE, the second fit uses 10 moments and 
applies the vector and axial- vector functions separately for the determination of as(m^) 
in combination with 5 parameters from OPE. Both fits yield nearly identical results for 
as(m^). The values are extrapolated to Mz using the RGE, the results for a;s(M|) are 
summarized in table |l[ The three different approaches describe the data equally well, the 
theoretical errors as well as the overall error is nearly the same, so from an experimental 
point of view there is no preferred calculation. However, the difference of the values 
measured applying different theoretical assumptions is about 4 % i.e. the difference is 
much larger than the the theoretical error determined for each method. The reason for this 
is most likely due to underestimation of the uncertainty due to missing higher order terms. 
For the ^^(m^) determination, the renormalization scale value = /x^/m^ has been fixed 
to x^ = 1, an uncertainty has been estimated due to scale variation in the range 0.4 < 

< 2, which is nearly the same range as in the DELPHI analysis of e"'"e~ event shapes. 
However, within the DELPHI analysis it turned out, that this scale variation range is 
sufficient only if one applies experimentally optimized scales, but yields inconsistent results 
for fixed scale values. An uncertainty due to the scheme dependence of the RGE coefficient 
Ps, which has been applied using the MS value, has been estimated due to variation 
between 0. < / 13^^^ < 2., however there is no theoretical reason, why its value should 
not be negative. This is indeed sometimes the case if optimized renormalization schemes 
are applied (se e.g. ||^). A similar analysis done by the ALEPH collaboration]^ has shown, 
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Figure 7: xFs as a function of x at the four lowest values. The curve shows a power 
law fit to the points with x < 0.1, which is used to calculate the integral in the shaded 
regions. 



that the use of the PMS scheme optimization leads to an reduced value of q;s(M|) and 
therefore reduces the discrepancy between FOPT and RCPT. Since the renormalization 
scheme optimization turns out to be of major importance in e^e~ annihilation, it would 
clearly be desirable to do similar studies in the analysis of hadronic r-decays. 

3.3 Oig Determination from the Gross-Llewellyn- Smith Sum 
Rule in Z/-N-DIS 

The Gross-Llewellyn-Smith (GLS) Sum Rule predicts the integral over the non-singlet 
structure function xF^lxjQ"^) measured in z/-N deep inelastic scattering. In the naive 
quark parton model, the value of this integral should be three. QCD corrections have 
been calculated in NNLO: 



1 



dx 



x 



xF,{x,Q')- = 3 l-^-a{nf){^y~b{nf){^y - — - (4) 



TT TT TT 



AHT 



The recent CCFR/NuTeV analysis |jTO[ covers an energy range of 1 GeV < < 
15 GeV^ with {Q^) ~ 3 GeV^ The dominant error comes from the contribution of higher 
twist terms AHT, which have been estimated to AHT = 0.15 ± 0.15 GeV^ The GLS 
integral is evaluated using the xF^ data separately in different bins. For very low x 
the data have been extrapolated using a power law fit (see fig. |^. Evolving the fit result 
for a,(3 GeV^) to M| yields = 0.114 ±:°°5 (^stat) ±;°°^ (syst) ± .005(t/ieo), which 

is consistent with other precise a.^ determinations in different energy ranges. 
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3.4 O-s Determination in Z/-N-DIS from high x scahng violations 



Another determination in I/-N-DIS within an updated analysis of the CCFR collab- 
oration |[Tl[] comes from a simultaneous fit of the F2 and xF^ structure functions. The 
QCD predictions are known in NLO. Improvements w.r.t. previous analyses are due to 
the new energy calibration and the higher energy and statistics of the experiment. The 
energy range of the data used for the as determination is 5 GeV^ < < 125 GeV^. The 
fit yields a,(M|) = 0.119 ± 0.002(exp) ± Omi{HT) ± 0.004(scaZe). The precision is 
better than from the NNLO analysis applying the GLS sum rule. 

The result can e.g. be compared with a result of similar precision from the analysis 
of the F2 structure function data from SLAC/BCDMS|T^. Within this analysis has 
been determined to as(M|) = 0.113 ± 0.003(exp) ± 0.004:{theo). Within both analyses, 
the renormalization scale value has been chosen to = 1. For an estimate of the 
scale uncertainty, has been varied within a large range of 0.1 < < 4.0. A similar 
variation has been done also for the factorization scale uncertainty. The range for the 
scale variation has been chosen in such a way, that the of the fit is not significantly 
increased, indicating that the choice = 1 for the central result is appropriate for the 
analysis of structure functions in DIS. 



4 Running of cXs 

The running of the strong coupling is among the most fundamental predictions of 
QCD. Due to the large energy range covered by LEPl/2, this prediction can now be 
tested from e"'"e~ data, measured within a single experiment. Whereas the perturbative 
predictions lead to an approximately logarithmic energy dependence of event shapes, the 
hadronization process causes an inverse power law behavior in energy and can therefore be 
disentangled from the perturbative part by studying the energy evolution of event shapes. 
Traditionally, hadronization corrections are calculated by phenomenological Monte Carlo 
(MC) models, whose predictions can now precisely been tested over a large energy range. 
Although MC-models yield a good description of the measured data within a large energy 
range (see e.g. [1^), their predictive power suffers from a large number of free parameters, 
which have to be tuned to the data (e.g. 10-15 main parameters to be tuned within the 
JETSET partonshower model). A novel way in understanding the hadronization process 



has been achieved with the Dokshitzer- Webber (DW) model|jT3[, which has recently been 
improved by results from 2-loop-calculations. Within this model, the power behaving 
contributions can be calculated, leaving only a single non-perturbative parameter to be 
determined from the data. 



4.1 Power Corrections and the Dokshitzer- Webber model 

Non-perturbative power corrections in the spirit of the DW-model arise from soft gluon 
radiation at energies of the order of the confinement scale. The leading power behavior 
is quantified by a non-perturbative parameter ao, defined by 
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ao(/ii 



1 



as{k)dk 



(5) 



which is expected to be approximately universal. Here, the true coupling as{k) is assumed 
to be infrared regular and can be understood as the sum of two terms 



(6) 



where the perturbative part a^^ is separated from the non-perturbative part a^^ by an 
infrared matching scale /x/ of the order of a few GeV. It is expected, that the factorial 
growing divergence of the fixed order perturbative expansion gets cancelled due to the 
merging with the non-perturbative counter-part, yielding a renormalon free theory. 



For the shape observables Y studied so far, the DW-model predicts the effect of the 
power corrections to be a simple shift of the perturbative distribution, i.e. for Y = 

1-T,C,Mh: 



da 



{Y - CyV) 



dY^ ' dY 

with an observable dependent factor cy and V proportional 1/Q: 



(7) 



V 



TT^ Q 

For the jet broadening observables Y 



— (Y) 
dY^ ' 



Ztt \ III (3q 
Bw, Bt the shift is predicted to be 



dY 



(Y 



Bo 



(9) 



with an additional non-perturbative parameter a'^ and a term V' entering via the log- 
enhanced term. Within the re-analysis of the JADE-data|]l5l it turned out that in order 
to describe the transition from the perturbative predictions to the observed spectra of the 
jet broadening observables, not only a shift, but also a squeeze of the partonic distribution 
is required. The original calculations for the jet broadening observables are considered to 



be erroneous, however the problem seems to be solved now[16 . 



4.2 Power Corrections to mean event shapes 

The earliest predictions from the DW-model have been made for the mean values of 
event shape distributions. In the context of the analysis of LEP2 data they have the 
advantage to make use of the full data statistics. Fig. | shows QCD fits in 0{a1) applying 



the DW-model to mean event shapes measured at various Ecm-, done by DELPHI [14 
and ALEPHU^. The fit results are listed in table |[ For the two observables studied, 
universality of ao is found on a 10 % level. Apart from a somewhat poor fit quality 
for the (1 — T) data used by the DELPHI collab., which can be explained by the poor 
quality of the low energy data, the DW-model fits yield a good description of the data, 
which is even better than for MC hadronization corrections. The as(M|) values obtained 
applying the DW model are in remarkable good agreement with as(M|) as determined 
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Figure 8: left side: Measured mean values of 1 — T and Mh from DELPHI high energy 
data and low energy data from various experiments. The solid lines represent the QCD fits 
in O(a^) applying the DW-model, the dotted lines show the perturbative contributions 
only, right side: (1 — T) distribution from an ALEPH analysis. The solid line represents 
a QCD fit in 0{a^) applying the DW-model the dashed line represents a fit applying MC 
hadronization corrections. 

by the DELPHI analysis of e"'"e~ event shapes at the Z^. There is however a fundamental 
difference between both analyses: Whereas the analysis of Z° data applies experimentally 
optimized scales in combination with MC hadronization corrections, the analysis applying 
the DW-model applies fixed scale values = 1, since a scale optimization is not feasible 
in this case. However, the fit quality of the DW-model fits indicate, that = 1 is an 
appropriate choice. The ALEPH result for the QCD fit to the (1— T) distribution applying 
MC hadronization corrections is also listed in table 0. Here, again = 1 has been applied. 
The fit quality is quite poor and as(M|) is much larger. The same observation has been 



Analysis 


Observable 


ao{2 GeV) 


«s(M|) 




DW-model (DELPHI) 


(1-T) 


0.494 ±0.009 


0.1176 ±0.0057 


43/22 


DW-model (DELPHI) 


{Mh) 


0.558 ±0.025 


0.1172 ±0.0037 


2/11 


DW-model (ALEPH) 


(1-T) 


0.497 ±0.009* 


0.1184 ±0.0019* 


28/30 


MC-corr. (ALEPH) 


(1-T) 




0.1307 ±0.0005* 


66/31 



Table 2: Results of QCD fits in O(a^) to the mean values of event shapes at different 
Ecm- Hadronization corrections have been applied either by the means of the DW-model 
or as predicted by MC. * Statistical errors only. 
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OPAL 183 GeV 
OI'AL 172 GeV 
ALEPH 172 GeV 
OPAL 161 GeV 
ALEPH 161 GeV 
OPAL 133 GeV 
ALEPH 133 GeV 
OPAL 91 GeV 
ALEPH 91 GeV 
DELPHI 91 GeV 
L3 91 GeV 
SLD91GeV 
JADE 44 GeV 
JADE 35 GeV 



Figure 9: left side: Measured distributions of Thrust, C-Parameter, wide jet broadening 
and the heavy jet mass from ALEPH data between 91.2 GeV < Ecm < 183 GeV. The 
hues represent a QCD fit applying the DW-model to shape distributions in matched 
0{a1) and NLLA precision, right side: QCD fits to the C-Parameter in matched 0{al) 
and NLLA precision applying DW power corrections to JADE and LEP data between 
35 GeV < ^^cm < 183 GeV. 



made within the DELPHI Z° analysis, and the result is just another demonstration, that 
fixed scale values are not qualified for an analysis applying MC corrections. 

4.3 Power Corrections to event shape distributions 

ALEPH ||T^ has studied event shape distributions of Thrust, C-Parameter, wide jet broad- 



ening and the heavy jet mass in matched 0{al) and NLLA precision. The central results 
are quoted as averages between two different matching schemes. They studied as well 
power corrections as hadronization corrections from MC. Fig. |^ shows the distributions 
in the energy range 91.2 GeV < Ecm < 183 GeV. For the fits applying the DW-model they 
reported a poor fit quality for the jet broadening distribution and that no determina- 
tion was possible. This can be understood due to the erroneous theoretical calculation 
mentioned before. For the heavy jet mass they got reasonable results applying the In R 
matching scheme but the fits were unstable under systematic variations of the analysis. 
The results for the Thrust and the C-Parameter fits are listed in table ^. 

As for the DW fits to mean event shapes, the quality of the fits to event shape dis- 
tributions is better if the DW-model is applied than for fits applying MC corrections. 
The agreement with DW-model fits to mean event shapes in 0{al) is good, in particu- 
lar for the combined fit to the Thrust and C-Parameter. It should be emphasized, that 
the observed agreement is different than the observation within the DELPHI analysis 
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A n 1 \r^\ ^ 
Xi.iicxi y oio 


O h cipyArpi hip 


(yn(2 CyeV) 

LtQl V.' ^ V 1 


n (Ml) 




T 


0.462 ± 0.060 


0.1193 ±0.0064 


DW-model 


c 


0.449 ± 0.082 


0.1130 ±0.0046 




C,T comb. 


0.451 ±0.061 


0.1168 ±0.0044 


MC-correction 


T 




0.1253 ±0.0061 


(LEP 1) 


C 




0.1212 ±0.0065 


MC-correction 


T 




0.1282 ±0.0054 


(LEP 2) 


C 




0.1233 ±0.0060 



Table 3: Results for the ALEPH fits to Thrust and C-Parameter. A combined fit to 
Thrust and C-Parameter has been made for the DW-model analysis. The fits applying 
MC hadronization corrections have been done separately for LEPl and LEP2 energies. 



of Z° data, where it turned out, that apart from a poor description of the data, the 
values from matched predictions were systematically higher than for 0{a1) predictions. 
The results from the DELPHI analysis suggested a mismatch between the quite different 
renormalization scale values required for NLLA and 0{a1) theory. As shown in the pre- 
vious section, a fixed renormalization scale value = 1 is appropriate if the DW-model 
predictions are applied, therefore there should be no mismatch if the matched predic- 
tions are applied in combination with the DW-model. The agreement between 0{a1) and 
matched results can then be interpreted in such a way that the contribution of the higher 
order logarithmic terms is quite small. In contrast the values from fits applying MC 
hadronization corrections are quite large, they are indeed larger than in any of the high 
precision analyses introduced before. This observation is basically the same than in |jl|], 
also the fit quality is worse than for the DW-model, but acceptable within this analysis. 



DW-model fits to event shape distributions have also been made within the re-analysis 
of the JADE data[|T3] between 35 GeV and 44 GeV in combination with various data from 
LEP experiments. Since the predictions for the jet broadening observables turned out to 
be erroneous, the results for the fits to Thrust and C-Parameter only are given in table 
^. So far, only the statistical errors have been evaluated. The results are in agreement 
with the results from the ALEPH collaboration. There seems to be a trend that the 
measured as values are even smaller than as from DW-model fits in 0{a1), contrary to 
as measurements applying MC corrections. 



Observable 


ao(2GeV) 


«s(M|) 


T 


0.501 ±0.009 


0.1136 ±0.0015 


C 


0.482 ± 0.008 


0.1128 ±0.0022 



Table 4: Results from QCD fits in matched 0{a^^ and NLLA precision applying power 
corrections predicted by the DW-model to the distributions of Thrust and C-parameter. 
The errors stated are statistical errors only. 
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Figure 10: Moments of the C-Parameter Distribution at different Ecm from an OPAL 
analysis. Superimposed is a fit of QCD including DW power corrections for the 

first moment and MC-corrections for the higher moments. The solid line shows a fit with 
experimentally optimized renormalization scales, the dashed line with = 1 kept fixed. 

4.4 Ois and its running from the higher moments of shape dis- 
tributions 

Even at LEP energies, perturbative QCD predictions yield a significant correction due to 
large non-perturbative, power suppressed corrections. This contributions can be reduced 
by the study of the higher moments n of event shape distributions: 



which has recently been performed within an OPAL analysis |jT8|. For shape distributions 
with power corrections proportional fJ^i/Q, the power corrections of the corresponding 
higher moments are expected to be suppressed by factors {ni/QY- The size of the power 
corrections could in principal be determined by applying the DW-model, in practice how- 
ever it turns out, that the corresponding non-perturbative parameter a^-i can only be 
constrained by the data for the first moment n = 1. Therefore DW power corrections 
have only been applied for the first moments and MC corrections for the second and third 
moment of thrust and C-Parameter. QCD fits in O(a^) have been done to the first three 
moments of event shape distributions simultaneously at various Ecm, Fig. [TO shows for 



example the moments of the C-Parameter. The results of the determination of as(M|) 
are listed in table ^. As for the analysis applying MC corrections introduced before, also 
OPAL finds only a poor description of the data, if the renormalization scale value is fixed 
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Observable 


Moments 


aJMl) 


/ Tlrlf 

A / '"Of 


a J Ml) y'^/n^f 






Xn fitted 




= 1, fixed 


C 


1-3 


0.116 ±0.010 


9.4/9 


0.141 ± 0.010 "large" 


1-T 


1-3 


0.114 dz 0.009 


5.8/9 


0.140 ± 0.009 "large" 


C 


1 


0.1164 


2.3/3 


0.1307 2.8/3 


C 


2 


0.1164 


2.7/3 


0.1537 3.1/3 


c 


3 


0.1164 


2.9/3 


0.1609 3.1/3 



Table 5: Fit results of QCD fits in Oip?^ to the moments of event shape distributions at 
energies between 91.2 GeV < E^m < 172 GeV. DW power corrections have been applied to 
the first moment only. Results on as(M|) are shown for simultaneous fits to the first three 
moments of 1-T and the C-Parameter applying experimentally optimized scales and for 
fits applying a constant renormalization scale value = /i/Q = 1. For the C-Parameter, 
also the as values for separated fits to the individual higher moments are given. For these 
fits, the uncertainties on as{M\) are not available so far. 



to = 1, whereas the description is perfect in the case of experimentally optimized 
scales. In the case of experimentally optimized scales, the fitted values of as(M|) are in 
good agreement with the results of the determination of as(M|) in C(a^) at the Z°|l[] 
as well as with the as(M|) determinations applying DW power corrections in 0[o?^ and 
matched 0[o?^ with NLLA. The largest contribution to the total uncertainty of as(M|) 
is due to the variation of the renormalization scale. The quite conservative estimate of this 
uncertainty yields Aas(M|) = 0.0074(sca/e) for both observables. For the fits applying 
a fixed scale value, the as(M|) values obtained are however much larger, and even under 
consideration of large scale uncertainties only poorly compatible with the precise as(M|) 
determinations introduced before. The differences of the results between experimentally 
optimized scales and fixed scale values are even more obvious, if one looks at the results 
for the fits to the individual moments of the shape distributions. They are listed for the 
C-Parameter as an example in table |^. Whereas in the case of experimentally optimized 
scales as(M|) is exactly identical for the separate fits to the first three moments, the 
scatter of as(M|) is about 20 % for the three fits. The largest value of as(M|) = 0.1609 
is obtained from a fit to the third moment, which differs from the current PDG-average 
by several a . A model with constant has been excluded within the OPAL analysis 
with a confidence level of at least 95 %. 



4.5 Running of ols from LEP data between 30 and 183 GeV 

Within the L3 analysis] TP] the running of has been studied including also e"^e~ anni- 
hilation data at reduced center of mass energies Ecm < M^. For this purpose, events with 
hard, isolated photons have been selected. This high energy photons are radiated either 
through initial state radiation or through quark bremsstrahlung, which takes place before 
the evolution of the hadronic shower. The QCD scale is assumed to be the center of mass 
energy of the recoiling hadronic system V^, which is related to the photon energy by 
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Figure 11: ag measurements from L3 data as a function of Ecm- The solid and dashed 
hues represent fits with an energy evolution as predicted by QCD and with constant a^, 
respectively. 




The selection provides determinations at reduced a/s^ from 30 to 86 GeV. ag has been 
determined from a QCD fit in matched 0(0^) and NLLA precision to the distributions 
of thrust, heavy jet mass, wide and total jet broadening. Hadronization corrections have 
been calculated by the means of MC models. Although the matched predictions turned 
out to be less reliable |l[] for the determination of an absolute value for as{M^), this 
should be no problem in terms of the running of the strong coupling, since the errors are 
fully correlated. Fig. shows the as(M|) values measured as a function of Ecm together 
with a fit to the QCD evolution equation. The fit yields a of 16.9 for 10 degrees of 
freedom, corresponding to a confidence level of 0.076, whereas a model with constant 
yields a of 91.4 corresponding to a confidence level of 2.9 x 10"^'^. 



5 Status of the Strong Coupling 

Measurements of the strong coupling are available from a large number of different reac- 
tions. Some problems arise, when the individual results are combined in order to calculate 
a global average. First, a global average of as(M|) contains a certain subjective element 
in the way the input data are selected. There are for example a large number of measure- 
ments in e"'"e~ annihilation at different center of mass energies which are however expected 
to be strongly correlated. The fact, that the exact correlation pattern between different 
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measurements is unknown suggests a pre-clustering of the input data in order to achieve a 
balanced mixture of measurements from different reaction types, which are then hopefully 
less correlated. Further problems arise due to the dominance of theoretical uncertainties 
within the determinations. Most experiments try to calculate uncertainties due to 
missing higher order contributions by the means of the variation of the renormalization 
scale, however, the range within the scale should be varied is quite arbitrary and different 
for each experiment. Therefore the resulting uncertainties on a<i(M|) are arbitrary to a 
large extent as well. 

Therefore, three different numbers will be given within the following considerations. 
The first number is a simple unweighted mean, which has within this context the advan- 
tage to ignore the doubtful scale ambiguity errors at all, however, different experimental 
uncertainties are ignored as well. The second number will be a simple weighted average, 
which does not account for the (unknown) correlations between the different mea- 
surements. Also an estimate for a correlated weighted average as introduced inpOf will 
be given. This method tries to estimate the covariance matrix by assuming a common 
correlation between the measurements, described by a single parameter pefr- The method 
applies, if < na/. Then, the measurements are assumed to be correlated and pes is 
adjusted until the expectation = ^df is satisfied. In general, this method yields a con- 
servative error estimate. The uncertainty of the average value might be quite large, if (e.g. 
theoretical) uncertainties are overestimated for some of the measurements included. Too 
small errors for the individual measurements yield in an underestimation of the correlation 
between the measurements. 



5.1 Status of cXs (early 1995) 

In order to illustrate the enormous progress achieved within the last three years, this 
overview is started with a summary of as(M|) measurement s from 1995. Fig. |12] 



shows a graphical overview of the different measurements. There were apparently two 
problems with the measurements shown. First the as measurements from Lattice Gauge 
Theory (LGT) ( as(M|) = 0.115±0.003) and from hadronic r-decays ( as(M|) = 0.125± 
0.004), which claimed both to be the most precise, yielded a large difference. They have 
been ignored in the global average, motivated by the fact that the LGT value has been 
unstable in the past and the r-decay value due to the controversial discussion about 
the validity of some specific theoretical assumptions. The second problem was, that the 
as(M|) values were clustered into two different groups of low and high energy {Q > 
25 GeV) measurements (with the exception of as(M|) from r-decays). The difference 
observed gave some input for speculations about new physics occurring at this energy 



value 


"s(M|) 




simple mean 


O.IU 


S5 ± 0.0072 




simple weighted average 


O.IU 


BO ±0.0022 


6.4 / 10 


correlated w. average 


O.IU 


BO ± 0.0045 





Table 6: A global average for as(M|) measurements representing the status in early 
1995, estimated using three different methods. 
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pp — > bB+jets 

e e [ev.shapes 34 GeV] 
e*e"[ev.shapes 58 GeV] 
pp — >W+jets 

Rz 
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0.04 0.06 0.08 0.1 0.12 0.14 



Figure 12: Summary of as{M'^) measurements at the beginning of 1995 [21]. The mea- 
surements are ordered according to the energy of the process. The shaded band represents 
the uncertainty of the correlated weighted average. 



scale. The overall scatter of the individual measurements is quite large (see table the 
uncertainty calculated from a simple weighted average clearly underestimates the true 
uncertainty, however the correlated weighted average seems to be appropriate here. 



5.2 Status of cts today 

In comparison with 1995 the situation now has drastically changed (see Fig. |TB| for 
a graphical overview on current as measurements). The difference between as(^|) 
from LGT and from hadronic r decays has been largely reduced. Since the q;s(M|) 
measurements from r decays revealed a somewhat larger uncertainty due to the differences 
from the models employed, a preliminary average of as(M|) = 0.119 ± 0.003 has been 
considered, corresponding to an average over the three different models introduced before. 



value 


«s(M|) 




simple mean 


0.1187 ±0.0022 




simple weighted average 


0.1189 ±0.0014 


1.84 / 10 


correlated w. average 


0.1189 ±0.0039 





Table 7: A global average for as(M|) measurements representing the status after the 
ICHEP'98 conference. The light shaded band represents the uncertainty of the correlated 
weighted average. 
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Figure 13: Summary of as(M|) measurements. Status after the ICHEP'98 conference. 



The current PDG average value for a,(M|) from LGT of as(M|) = 0.117 ± 0.003|22[ is 
in good agreement with the value from r decays and there is no longer a reason to exclude 
them from the global average. Furthermore, the two clusters of values observed earlier 
completely disappeared. The global average for as{M^) (see table |^ ) is nearly the same 
than in 1995, however the scatter between the individual measurements is largely reduced. 
The is 1-84 for 10 degrees of freedom, indicating that the (theoretical) uncertainties 
might be overestimated. However, the procedure for calculating a correlated weighted 
average interprets the small entirely in terms of correlations between the observables, 
therefore yielding a quite conservative error estimate. The largest deviation from the 
global average comes from the event shape measurements in e"'"e~ annihilation. Here, the 
results from matched 0{a1) and NLLA predictions in combination with MC corrections 
have been used for the global average, since this is the standard method nowadays and a 
large number of measurements have been made. But as we have seen within the previous 
discussion, there are serious arguments, why this predictions yield less precise results 
than predictions in Ola^) in combination with experimentally optimized scales. It is 
quite instructive to study the changes on the global average value, when the results from 
matched predictions are replaced. This will be done in the next section. 

5.3 looking into the future ... 



In a first step, the global results from data have been replaced by the single 
result of as(M|) = 0.117 ± 0.003 from the DELPHI collaboration 0, obtained from 
0{al) predictions in combination with experimentally optimized scales. Secondly, it has 
been assumed, that the DW-model gets estabhshed and the results from e+e~ data with 
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Figure 14: Another summary of as(M|) measurements. Here the results from matched 
C(q!^) and NLLA predictions have been replaced with results from 0{a^) QCD applying 
experimentally optimized scales, and with results from using DW-model corrections. The 
light shaded band represents the uncertainty of the correlated weighted average. Addi- 
tionally the uncertainty of the simple weighted average is indicated by the dark shaded 
band. 



Ecm 7^ Mz have been replaced by an average value of current results from DW-model 
predictions. This average has been calculated assuming fully correlated errors which 
yields as(M|) = 0.117 ± 0.005. (Here, the results of |§ have been ignored, since only 
statistical errors are given.) See Fig. |l^ for a graphical view and table |^ for the average 
value. The modified result is quite impressive. The global average is about 1 % smaller 
than before, and the consistency of the measurements gets further improved. The scatter 
of the individual measurements is now only ±0.0012 and the uncertainty of Aa^l^f) 
= ±0.0038 seems now really be too pessimistic. The uncertainty obtained from a simple 



weighted average is about 1 % and indicated in Fig. |1J for illustration reasons. Clearly, 



no correlations between the measurements are considered here. Also the new results 



value 


«s(M|) 




simple mean 


0.11 


B5± 0.0072 




simple weighted average 


0.11 


BO ±0.0022 


6.4 / 10 


correlated w. average 


0.11 


BO ± 0.0045 





Table 8: A modified global average for a;s(M|) where the results from matched predictions 
have been replaced. 
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introduced here have to be estabhshed. However, if one considers the progress achieved 
within the last three years, a 1% error on as(M|) seems to be a reahstic perspective for 
the foreseeable future. 



6 Summary 

An enormous progress has been achieved on the determination of a<j(M|) and its running 
with the analyses presented at this years summer conferences. The importance of adjust- 
ing the renormalization scale has been demonstrated with the analysis of high statistics 
and high precision TP data using angular dependent shape observables. The observation 
is confirmed by the analysis of higher moments of event shapes from LEP high energy 
data. Comparison of the data with predictions from matched 0{a1) with NLLA preci- 
sion revealed an so far unreported problem, presumably arising due to a mismatch of the 
renormalization scales. The DW-model is a novel way in understanding non-perturbative 
power corrections to event shape distributions. First results are impressive. Applying the 
DW-model, QCD predictions in 0{a1) precision and matched 0{a1) and NLLA precision 
yield similar as{M'^) values, which are also in good agreement with other results from 
precise as{M'^) measurements. If hadronization corrections from MC models are applied 
instead, a larger deviation of as(M|) is observed for matched predictions, however still 
compatible with the global a<((M|) average value. The running of as is clearly established 
from LEP e"^e~ data only. All recent precise as(M|) measurements agree very well with 
each other, a global average of 



a,(M|) = 0.1189 ±0.0039 

has been determined using a correlated weighted average. Replacing the results from 
matched 0{a1) and NLLA predictions with the new results obtained from applying ex- 
perimentally optimized scales and from DW model predictions further improves the con- 
sistency of the global as(M|) measurements. A roughly 1 % error on as(M|) seems to 
be a realistic perspective for the foreseeable future. 
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